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Qh| its impact on — > K*^K^ decays and their CP-conjugate processes. We show that 
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resonances at a r-charm factory, it is even possible to determine or constrain x, y and 
5 from the time-independent measurements of coherent (D^D^) {K*^K^){K*^K^) 
' decays. If the CP-violating phase of D^-D^ mixing is significant in a scenario beyond the 

■ standard model, it can also be extracted from the K*^K^ events. 
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1 1 1 The BaBar [1] and Belle [2] experiments have recently provided us with some noteworthy 
evidence for D^-D^ mixing, a quantum phenomenon similar to K^-K^, B^-B^ or B^-B^ mixing. 
Both experiments indicate a non-vanishing width difference between the mass eigenstates 
and D2 of D° and mesons, 

y' cos = (0.97 ± 0.44 ± 0.31) x 10"^ , 

ycp = (1.31 ± 0.32 ± 0.25) X 10-2 , (1) 

where the values of y' cos (f) and y^p are extracted from the decay modes K^tt~ 
versus K+n- [1] and ^ ^ K+K' and n+n- versus D'^ K+R- and tt+tt- [2], 

respectively. By hnking y' and y^jp to the D^-D^ mixing parameters x = (M2 — M-^)/T and 
y = {T^- ri)/(2r) with r = (F^ + 12)72 and T^ ^ being the width of D^ ^, Nir has pointed 
out that \y\ ~ 0.01, \x\ < \y\ and small or vanishingly small CP violation are expected in the 
D^-D^ mixing system within the standard model [3]. Some other authors have also discussed 
possible implications of the BaBar and Belle measurements of D^-D^ mixing, either within or 
beyond the standard model [4]-[9]. 

Unfortunately, current theoretical calculations of D^-D^ mixing involve large uncertainties 
because of the dominance of long-distance contributions [10]. In the standard model, the values 
of X and y are expected to be a second-order effect of the S'f/(3) flavor symmetry breaking [11]: 
x,y ^ sin^^^^ x [5't/(3) breaking]^, where 6q ^ 13° denotes the Cabibbo angle. A very reliable 
prediction for the size of SU (3) breaking has been lacking, although many attempts have been 
made [10, 12]. Hence these two D^-D^ mixing parameters might be only of limited use in 
testing the standard model and searching for new physics. From an experimental point of view, 
however, it is always desirable to measure or constrain x and y as accurately as possible. 

Motivated by the afore-mentioned positive results from the BaBar and Belle experiments, 
here we aim to investigate the impact of D^-D^ mixing on — > K*^K^ decays and their 
CP-conjugate processes. Because K*^K'' (or K*''K'^) is not a CP eigenstate, the amplitudes 
of D° and decays into K*'^K'' (or K*~K'^) may have a significant strong phase difference 
5. In contrast, vs — > K^K~ decays do not involve such a strong phase difference. We 
show that both the D^-D^ mixing parameters {x and y) and the strong phase difference (5) can 
be determined or constrained from the time-dependent measurements of D*^ vs K*^K^ 
decays. On the ■?/'(3770) and ^(4140) resonances at a r-charm factory (e.g., BEPC-II [13]), we 
find that it is even possible to determine or constrain y and 5 from the time-independent 
measurements of coherent {D^D^) {K*^K^){K*^K'^) events. If the CP- violating phase of 
D^-D^ mixing is significant in a scenario beyond the standard model, it can also be extracted 
from the decay modes under discussion. 

The remaining part of this paper is organized as follows. Section 2 is devoted to the effects 
of D^-D^ mixing and CP violation in the time-dependent vs — > K*'^K^ decays. The 
coherent (D^D^) (X*±XT)(X*±XT) decays on the ^(3770) and ^'(4140) resonances are 
discussed in section 3, where we focus our interest on possible signals of D^-D^ mixing and CP 
violation in the time-independent events. An isospin analysis of the final state interactions in 
D — > KK* modes is done in section 4. Finally, we summarize our main results in section 5. 
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In the standard model K*^K^ transitions can occur through both tree-level and 

loop- induced (penguin) quark diagrams. The former is essentially CP-conserving (proportional 
to V^gV*g), while the latter is negligibly small (suppressed by m^/M,^ for q = d,s,b) [3]. 
Hence the four amplitudes = {K*+K'\n\D^), Aj^,+j^^ = {K*+ K-\n\D°) , Aj^,.j^+ = 



{K*-K+\n\D°) and = {K*- K+\'H\D^) have the relations Aj^,.j^+ 

Aj^,_j^+ = Aj^,+j^_ as a good approximation. We define 



Aj^*+j^— and 



^K*-K+ 
^K*-K+ 



pe 



iS 



(2) 



where p > and 5 is the strong phase difference. On the other hand, two neutral D-meson 
mass eigenstates can be written as 



(3) 

where p and q satisfy the normalization condition \p\'^ -f |gp = 1. The phase oi q/p is (f) = 
&Tg{q/p) = aTg[{V*gV^g) I ( V^csKTs)] ~ within the standard model [14] , but it might be significant 
if a kind of new physics contributes to the box diagram of D^-D^ mixing [3, 4, 5, 6, 12, 15]. 
Allowing for both \q/p\ ^ 1 and 7^ 0, we may use the following rephasing-invariant quantities 
to express the decay rates of D° K*'^K'^ and L>° K*'^K'^: 



q A 



K*+K- 



K'+K- 



K*-K+ 



P ^K*+K- 

q A 



P 



K*-K+ 



(4) 



Since the naive factorization approximation yields p ~ it is quite natural to expect that 



lA 



K*+K- 



iV-i^+l ~ C'(l) holds. 
First, let us look at the time-dependent decay rates of vs Z)° K*^K^. Now that the 
D^-D^ mixing parameters x and y are both small, we may just keep the terms of 0{x) and 
0{y) in our calculations. Using the generic formulas given in Ref. [16] ^, we explicitly obtain 



r[D''{t)^ K*+K-] <x 
r[D''(t)^ K*-K+] (X 



\^K*+K- 
\^K*-K+ 



'^e ^*[1 + {yReXf^*+f(- - xlmX 



2e-^*[l + {yRe\K*-K+ - xImA^.-^+)n] ; 



(5) 



and 

T[D\t) 
T[D\t) 



K*-K+] oc \A 



K*-K+\ 



e [\^K*-K+\ + {yReXK*-K+ + x^^^K*-K+)^t\ 



K*^K'] oc \A 



K*+K- 



■i-Note that y = — y^, where y^, is the D^-D^ mixing parameter defined in Ref. [16]. 



(6) 
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where we have required t <1/V for the proper time t. Taking account of Eq. (4) and defining 
the effective D^-D^ mixing parameters 



X COS 5 ± y sin 5 , 
y cos 5 ± a; sin 5 , 



(7) 



we simphfy Eqs. (5) and (6) to 

TlD'^it) ^ K*+K-] oc \AK.+K-?e- 

r[L)°(i) ^ X*-/s:+] oc \AK*-K+?e- 

and 





1+P 


1 




P 


rt 


1+P 


p 







i2„-rt 



T[D''{t) K*-K+] oc \AK*-K+\-e 







p' + p 








p' + p 


p 







[y'_ cos — sin 0^ Ft 
(^^/^ cos + x'j^ sin j rt 

{il'^ cos — sin 0^ 
(^1/^ cos + sin 0^ Tt 



(8) 



(9) 



Once these four decay rates are measured, it will be possible to determine p and constrain 
the magnitudes of both D'^-D^ mixing and CP violation. Note that the deviation of \p/q\ (or 
\q/p) from unity, which can also be determined or constrained from other neutral D-meson 
decays, signifies CP violation in D^-D^ mixing. This effect is conveniently described by a small 
parameter A up to the correction of O(A^); i.e., p/g = 1 + A and q/p = 1 — A. Given (p ^ Q 
in the standard model, useful information on y'j^ and y'_ is achievable from the time- dependent 
measurements of vs Z)° K*^K^ transitions. A clear difference between y'j^ and y'_ will 
imply that both x and 6 are not very small. These points have also been observed in Ref. [15]. 

We remark that the K*^K^ events of neutral D-meson decays are important, since they can 
be complementary to the X^tt^ and K^K~ (or 7r''"7r~) events for the experimental searches for 
both D^-D^ mixing and CP violation. A similar idea, which makes use of the D*'^D^ events 
of neutral S-meson decays to extract the CP- violating phase /3 and test the factorization 
hypothesis [18], has actually been adopted by the Belle [19] and BaBar [20] Collaborations in 
their experiments at the KEK and SLAC B factories. 



3 Now we turn to the possibility of measuring coherent {D^D^)q — >■ [K*^K^){K*^K^) 
decays on the ■?/'(3770) resonance with C = —1 and (or) on the ■?/'(4140) resonance with C = +1, 
where C denotes the charge-conjugation parity of the and pair. Both time-dependent and 
time-integrated rates of a general {D^D^)^ j\f2 decay mode, together with their approximate 
expressions up to the accuracy of 0{x'^) and C(|/^), have been formulated in Ref. [16] without 
special assumptions. Here we focus our interest on the time-independent measurements of those 
[K*^ K^)[K*^ K^) events from coherent {D^D'^)q decays at a high-luminosity r-charm factory 
(e.g., BEPC-II [13]). Let us define rj+ = r(K*+K-,K*+K-)c, Tj" = r{K*+K-,K*-K+)c, 
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= T{K*-K+, K*+K-)c and T^" = r{K*-K+, K*-K+)c for four joint decay rates. With 
the help of Ref. [16] ^, we exphcitly have 



r++ oc 2|A^ 



4 


p 


{ 




Q 






K- 








4 


1 


\ 




p 






K+ 



(2 + C) r 1 + C\\.+K- + (1 + Cy 



K*+K- 



ReA^.+^- +x(l- |A^.+ii--|^) ImA^.+ii-- } , 



(2 + C) r 1 + + (1 + C)' 



Ai<-* 



K*-K+\ 



+y('^ + \^K*-K+\'^)'^^^K*-K++x(l-\XK*-K+\^)'^T^^K*-K+]} '■, (10) 



and = with 



r+- oc 2\Aj,.+j,-\U{2 + C)r 



+ 



1 + CXf^t+T^-Xf^t.- 



^K*+K-^K'—K+ 
*+K- 



*+K- *-K+\ 

+ (1 + C) y [(l + \\K.+j^-f) ReXK.-K+ + (l + |A;^*-;^+|^) ReA;^, 
-(1 + C)a;[(l- |A;^*+^-|2)lmA;^.-;^+ + (l- |A;^.-;^+|2)lmA^.+^-]} , (11) 



where r = (x^ + y^)/2 is essentially the ratio of wrong-sign to right-sign events of semileptonic 
and decays [16, 17]. When Eq. (4) is taken into account, Eqs. (10) and (11) can be 
simplified to 



r++ oc 2\Aj,.+j,.\''{{2 + C)t 



1 


2 


p 





+ (i + c)V 



r^- oc 2\A^.-^^\H{2 + C)r 



+ (1 + C)V 



+ 2Cp^cos((5 + 0) 
{^'j^ cos -I- x'_ sin (j)^ + - {y'__ cos — x'^ sin 0^ | , 

+ 2Cp2 cos (5 - 0) + 
(^y^ cos (f) — x'_ sin (j)^ + p"^ - {y'_ cos -|- x'_^ sin 0^ | ; 





2" 


p 




Q 





(12) 



and 

r 



c 



oc 2|^;^,+;^_n (2 + C)rp2 



+ 2Ccos(20) + 



+ 



1 + 2Cp2 cos (25) + 



+ {1 + C)p 
+ {l + C)p 



p 









[y'_ cos <p + x'j^ sin j -|- p^ cos -|- x'_ sin j 



cos — X, sm 



+ cos — x'_ sin j 



(13) 



Note that the terms proportional to r in r^"*" are only important when C = — 1 is taken. As 
K*^K~ and K*~ are not the CP eigenstates, both p 7^ 1 and 5 7^ are expected to hold. 



^Note again that y = —yj^, where yjj is the D°-D° mixing parameter defined in Ref. [16]. 
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It is therefore reasonable to neglect the term proportional to r in Eq. (13) even for the C = —1 
case. We stress that these formulas will be very useful to analyze the experimental data on 
coherent {D^D^)^ — * {K*'^K'^){K*'^K'^) decays at a r-charm factory. 
For C = — 1 on the ■0(3770) resonance, we obtain 



rl+ _ ^1 -2p2cos(5 + (/)) +p' 



ri" 1 - 2p2 cos (25) + p4 

n- _ 1 - 2p2 COS {5 -(()) + 
fl^ ~ ^ 1 - 2p2 cos (25) + p4 



(14) 



where we have used the approximation \p/q\ ~ \q/p\ ~ 1 and neglected the term proportional 
to r in rl;~. One can clearly see that these two ratios signify D^-D^ mixing (i.e., r 7^ 0). The 
difference between Fl^/rl" and ri~/rl~ measures the CP-violating effect in D^-D^ mixing 
(A ^ 0) and that from the interference between decay and mixing {(f) ^ Q): 



r++ _ p-- 4r (1 — p'') A + p^ sin 5 sin 1 



1 - 2p2 COS {25) + p4 



(15) 



where the notations \p/q\ = 1+A and \q/p\ = 1— A have been taken into account. The smallness 
of r (i.e., r ~ 10^^), however, might more or less obstruct the observation of Tt.'^/T^'^ and 
ri^/r+^ at present. But we hope that the high-luminosity r-charm factory may finally realize 
the desired measurements in the near future. 

For C = +1 on the 0(4140) resonance, one may simply neglect the terms proportional to r 
in Eqs. (12) and (13). Up to small corrections of 0{x'^) and C(2/±), the relationship 

r++ 4^2 

r+- ~ r+- ^ 1 + 2p2 cos (25) + p4 ^ ' 

holds approximately. Once the ratios r;}:+/r:}:~ and T';^" /T'^^" are measured, they will impose 
a strong constraint on p and S. On the other hand, the difference between r|;+/r;|;~ and 
r;;j:~/r:[;~ is a clear signal of CP violation: 



r++ 



A cos (y'j^ - p^y'S) + sin (y'_ - p^x'^ 



Vf 1 + 2p2 cos {25) + p^ 



(17) 



Comparing between (rl+ — ri~)/rl^ and (F^"*" — F^~)/F^~, we find that the latter is less 
suppressed by the smallness of x and y. Hence it seems more promising to measure CP violation 
in the decays of correlated and mesons into {K*"^ K^){K*^ K"^) states on the 0^(4140) 
resonance. 



Finally let us make some comments on the final-state interactions in — > K*^K^ 
transitions. A model-independent approach is to do the isospin analysis of — > K*^K', 
D° K*^K^ and D+ K*+K° decays or K*-K+, K*^K^ and D+ K*^K+ 

decays, in which each final state contains 7 = and (or) I — 1 isospin configurations. For 
simphcity, we denote the amplitudes of D° K*+K-, D° ^ K*^K^ and D+ K*+K^ as 



6 



*^K*+K- 



.0/^0 Slid Aj^^ 



respectively. They can be expressed in terms of two independent 



isospin amplitudes and ^4^^ as follows [16]: 



2(^+^0) , 

I (A - A) , 
A- 



(18) 



The branching ratios of these three decays are B^,+j^^ oc |A^,+^_prQ, -B^.o^^o oc |^^*o;^opTo 
and B^^+j^o oc |A^,+^o|V+, where Tq = (410.1 ± 1.5) x 10"^^ s and t+ = (1040 ± 7) x 10"^^ s 
are the hfetimes of and D'^ mesons [21], respectively. Defining ^IqMi — ^^^'^^ 

1 + ze' 





2 




^K*+K- 


1 










4 




2 


'[± 




1 






T"0 




4 



1 - ze'"^ 



(19) 



Then the isospin parameters z and (p can be determined: 



B 



K*+K 



-+B 



To 



B 



- 1 



1/2 



if — arccos 



B 



R'+R- 



B 



R*ORO 



zB, 



(20) 

'0 '^J^K'+R" I 

Of course, 7^ implies the existence of final-state interactions. One may follow a similar 
procedure to do the isospin analysis of K*^K^, K*^K^ and K*^K+ 

decays. The amplitudes of these three transitions are essentially identical to those of — > 
K*'^K~, — > K*'^K^ and D~ — > K*^K~ transitions, since their tree- level quark diagrams 
are CP-conserving and the penguin diagrams are negligibly small in the standard model. The 
corresponding isospin parameters z and (p can be extracted from the branching ratios 
^r*oro and -B^«ox+ the CP-conserving case. It is in general difficult to link the isospin phase 
differences (p and (p to the strong phase difference S defined in Eq. (2), unless some assumptions 
are made in a specific model of hadronic matrix elements. Nevertheless, it is reasonable to argue 
that significant (p and ip must hint at significant S for K*^K^ events. 

For the purpose of illustration, we do a numerical analysis of the isospin parameters by 
using the present experimental data [21], 

^.+^0 = (3.2±1.4) X 10-^ 
= (3.7±0.8) X 10-^ 



B 



Br*+r- 
Br *0K0 



< 1.6 X 10 



-3 . 



(21) 



and 



Bj^,OK+ = (3.02 ± 0.35) X 10" 
Bj^.-j^+ = (2.0 ± 1.1) X 10 



-3 



B 



R*ORO 



< 8 X 10 



-4 



(22) 
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Since the magnitudes of i?^,o^o and 5^,0^0 have not been fixed, our analysis can only provide 
some limited information on {z,Lp) and {z,(p). The numerical results are shown in Fig. 1 and 
Fig. 2. Some comments are in order: 

1. It is straightforward to sec that the possibility of ip = and (or) <^ = is almost excluded 
by current experimental data. The most favorable values of (/? and (p are around 50°, 
implying the presence of significant final-state interactions. Indeed, (f can be as large 
as 65°, and (/? can be even larger than 80°. The strong phase difference 6 is therefore 
expected to be significant in vs — > K*^K^ transitions. 

2. The constraints on B j^^q / B and -B^*o/fo/-B^*o^+ allow us to extract the lower 
and (or) upper bounds of -B^*o;^o and S^^o^fo- We find 3.2 x lO^'^ < B^^a^^ < 1.6 x 10~^ 
and 5^,0^0 < 8 X 10^^. The former is interesting and can be tested in the upcoming 
experiments, but the latter is trivial. More accurate data will reduce the uncertainties in 
our isospin analysis. 

3. The allowed ranges of z and z do not have much overlap. In particular, z > 0.8 > 2; is 
roughly true. This observation implies that K*^K^ and K*^K^ decays 
might involve quite different final-state interactions, from which significant 5 is naturally 
anticipated. 

It is worth mentioning that an isospin analysis of — > K~^K~, — > K^K^ and D'^ — > K^K^ 
decays [16], whose branching ratios have all been measured, also indicates the existence of strong 
final-state interactions. As K^K~ is a CP eigenstate, however, the ratio of {K^K~\H\D^) to 
{K^ K~\1-L\D^) does not involve a significant strong phase difference in the absence of direct 
CP violation [3]. 



In summary, we have investigated D^-D^ mixing and CP violation in — > K*^K^ 
decays and their CP-conjugate processes, whose final states may have a significant strong phase 
difference. We have shown that both the D^-D^ mixing parameters {x and y) and the strong 
phase difference {5) can be determined or constrained from the time-dependent measurements 
of vs — > K*^K^ decays. For a high-luminosity r-charm factory running on the ■0(3770) 
and ■0(4140) resonances, we find that it is even possible to determine or constrain x, y and 5 
from the time-independent measurements of coherent [D^D^) {K*'^K^)[K*'^K^) events. If 
the CP- violating phase of D^-D^ mixing is significant in a scenario beyond the standard model, 
it can also be extracted from the decay modes under discussion. 

We strongly recommend the experimentalists to pay some special attention to the K*^K^ 
events of neutral D-meson decays, because they are complementary to the K^n^ and K^K~ 
(or 7r"'"7r~) events for the study of both D^-D^ mixing and CP violation. We expect that these 
interesting channels and possible new physics in them can well be explored at BEPC-II and 
other charm-physics experiments in the near future. 

One of us (Z.Z.X.) would like to thank Phil Chan and the National University of Singapore 
for warm hospitality, where this paper was written. He is also grateful to H.B. Li for some 
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Figure 1: Numerical illustration of the allowed ranges of (/? and -B^«oi^o/-6j^«+^o- 
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